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Abstract 

The behaviors of one-dimensional quantum random walks are strik- 
ingly different from those of classical ones. However, when decoherence 
is involved, the limiting distributions take on many classical features 
over time. In this paper, we study the decoherence on both posi- 
tion and "coin" spaces of the particle. We propose a new analytical 
approach to investigate these phenomena and obtain the generating 
functions which encode all the features of these walks. Specifically, 
from these generating functions, we find exact analytic expressions 
of several moments for the time and noise dependence of position. 
Moreover, the limiting position distributions of decoherent quantum 
random walks are shown to be Gaussian in an analytical manner. 
These results explicitly describe the relationship between the system 
and the level of decoherence. 

©2007 United States Copyright Office 



1 Introduction 

Quantum random walks recently gained great interest from physicists, com- 
puter scientists and mathematicians. The interest was sparked by their im- 
portant roles in developing highly efficient quantum algorithms. For instance, 
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Grover's search algorithm pQ has time cost 0{-\fN), in contrast to the or- 
dinary search algorithm which has a cost of 0{N). This quantum search 
algorithm was proved to be closely related to the behaviors of quantum ran- 
dom walks in |[2j and |3] . As another example, Shor's algorithm also improved 
the speed of factorization dramatically The high efficiency of quantum 
algorithms is discussed in [3], [B] and [7]. Experimental implementations of 
the algorithms are discussed in [S] and [5]. 

Besides their important applications, quantum random walks are very 
attractive due to their dramatic non-classical behavior. After quantum ran- 
dom walks were defined in [10], many articles ([H], [12], [13], [E], [I5], etc.) 
studying the distribution of quantum random walks were published. It is 
known that the observed non-classical behavior is due to quantum coherence 
[16j. One of the most shocking differences fTT] is that as time t grows, the 
variances of quantum random walks are 0{t^) while the variances of classical 
random walks are 0{t). Various limit theorems of quantum random walks 
are established ([12], [IS], [IT], [H]). An excellent reference can be found in 
[19]. 

One of the most important issues surrounding the use of quantum random 
walks is that they are very sensitive to inevitable decoherence, which could 
be caused by many reasons, such as interactions with the environment and 
system imperfections ([20], [H], [22] )• The effect of decoherence is very 
important for the application of quantum algorithms, as discussed in [23], 
and [21]. For the one-dimensional case, in the model in [20], decoherence is 
introduced by measurements on the particle's chirality. Long-term first and 
second moments of the walk were obtained and numerical results showed that 
the distributions look like classical normal distributions. Similar results are 
found in other models ([IS], [23], [25], [2S], [2Z], [2B], [22]). If we denote by Xt 
the position of the particle obtained by measurement of the wave functions 
of the decoherent quantum random walk, all of above papers mentioned the 
fact that the variance of the simulated Xt grows linearly in time for large t 
when the quantum random walk is subject to decoherence. 

These results stimulated us to prove that the one-dimensional decoherent 
quantum random walk, converges to a normal distribution. Our work 
focuses on the one-dimensional discrete-time Hadamard walk with measure- 
ments taken on both position and chirality at each time step. This kind of 
decoherence is studied numerically in [8] , [23] , [25] and [26] but we will study 
it fully analytically. 

We shall see that when the particle is not measured, then the system 
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particle is measured subject to a small probability, will converge to a 
normal distribution. In the limit, when the particle is measured at each step, 
then the system becomes purely classical and ^ is asymptotically standard 
normal. 

In the next section, we introduce the mathematical setup of decoherent 
quantum random walks. We then provide our methodology of generating 
functions and the decoherence equation. We next list our results and discuss 
the interesting phenomena that occur when p is small. Finally, wc summarize 
and discuss our work. Mathematical proofs are given in the appendix. 

2 Decoherent Quantum Random Walks 

2.1 Pure Quantum Random Walk system 

We start with a brief description of the one- dimensional pure quantum ran- 
dom walk system. In the classical random walks, the particle moves to the 
right or left depending on the result of a coin toss. However, in the quantum 

random walks, the particle has its chirality {right, left} as another degree 
of freedom. At each step, a unitary transformation is applied to the chirality 
state of the particle and the particle moves according to its new chirality 
state. 

We denote the position space of the particle by Hp, the complex Hilbert 
space spanned by the orthonormal basis {\x),x G Z}, where Z is the set of 
integers. We also denote the coin space by He as the complex Hilbert space 
spanned by the orthonormal basis {|/),/ = 1,2} where 1 stands for "moving 
right" and 2 stands for "moving left." The state space H of the particle is 
defined as 



A vector e H with L^-norm 1 is called a state and tells us the dis- 
tribution of the particle's position and chirality upon measurements. The 
basis of H are denoted by {\x,l) = \x) ^ \l) : x ^ Z},/ = 1,2}. Now we 
introduce the evolution operator which drives the particle. The shift operator 
S:H^ His defined by 



H ^ Hp® He. 



(1) 
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The coin operator C : ^ He can be any unitary operator and is an 
analogue to the coin flip in the classical walk. The evolution operator U : 
H ^ H is defined by 

U = S{Ip0C), (3) 

where Ip is the identity in the position space. 

Now let |'0o) & H he the initial state and let |^t) = U^li^o)- The sequence 
(l^t)}^ is called a one- dimensional quantum random walk. 

The most famous and best-studied example of quantum random walks is 
the Hadamard walk, in which the coin operator is the 2x2 Hadamard matrix 

The quantum random walk associated with H2 is called a one- dimensional 
Hadamard walk. 

The probability of a particle at state l-^) to be found at state \r]) is defined 
by the norm squared of the inner product of and \r]), K'i/'lTy)^. 

In particular, the probability of the quantum random walk, starting from 
the position x = 0, with the coin in state m, to be found at x with coin state 
n is 

Wm,ni^,t) = \{x,n\U'\0,m)\^. (5) 



2.2 Decoherence 

We focus on decoherence caused by measurements that is on both position 
and coin of the particle. A set of operators {Ai,i e A} is called a measure- 
ment if 

J2a:a = i, (6) 

where A is some index set and A* is the adjoint operator of A, i.e., the 
complex conjugate of transposed matrix of A. 

In this work, we consider the measurements in a similar manner as in 
|20j . Let p be a real number in [0, 1], to denote the probability of the random 
walk being measured at each step. We let Ac : H ^ H be s.t. Ac = 
^/l — pi. Hence, Ac is the coherence projection. We also let A^^n . H ^ H 
be s.t. A^^ri = \fv\^i n) {x, n\ to be the decoherence projection to the subspace 
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spanned by \x,n). Under this setup, the index set A is A = {c} \J{{x,n) : 
X G Z,n = 1,2}. 

Let be a state in H. Then the position of the particle, Xf over 
Z obtained by taking measurement of the wave function of the decoherent 
quantum random walk starting from 1-0), at time t has the following distri- 
bution: 



In other words, the walk starts at \ip), then we apply the evolution oper- 
ator U and try to measure it. The process repeats until the t-th step is fin- 
ished. We then consider the position distribution of the particle. We call each 
(ji, j2, . . . , jt, (x,n)) a path. We also call {x,n\{AjJJ){Aj^_JJ) . . . {AjJJ)\il)) 
an amplitude function of the particle associated with the path. Many paths 
yield amplitude due to the decoherence projections, the A^j's. However, 
the summation in ([7]) over all paths {ji,j2, ■ ■ -jt, {x, n)) gives the probability 
of observing the particle at position state \x) at time t. 

At each step of a path, the walk is either not measured with probability 
q = 1 —p or is measured at n) with probability p. So when p = 0, the walk 
is not measured and the system is the same as the pure quantum random 
walk previously defined. When p = 1, the particle is interfered with at each 
step, hence the quantum behavior essentially disappears and the system is 
exactly classical. 

We work on the decoherent Hadamard walk starting from position 0. We 

use 



to denote the probability that, at time t, of a particle in the decoherent quan- 
tum random walk starting from 1-0) to be found at state In particular, 
we denote the probability that at time t, the particle starting at |0,m) can 
be found at \x, n) by 



P{Xf = x) 



(7) 







(9) 
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Since we are interested in the limiting distribution of the walk, we focus 
on the Fourier transform of the above probabilities, 

i^m,n(^,t) = 5^Pn^,n(x,t)e*'^ (10) 

X 

We consider two types of walks. We first consider the walk starting at the 
state |0o) = 1) + ^75|0'2). We call this walk "symmetric" and denote 
it by Xt. Note that the characteristic function of Xt is 

%'^) = ^EE^-.-(^'^)- (11) 

m n 

From the above equation, we can see that its characteristic function is ob- 
tained by taking the average of those with initial chirality state m. Further- 
more, in [11], it is shown that the pure quantum random walk starting with 
|0o) has a symmetric position distribution. These are the reasons why we 
call it "symmetric." 

We also consider the walk that starts at |0, 1), Xt, i.e, the walk starting 
at with chirality "right." In this case, the characteristic function of Xt is 

P{k,t) = J2PiAk,t). (12) 

n 

Our goals are to show that as t ^ oo, 

P(A;,t) ^6^5-'=', (13) 

for some positive number v in the symmetric walk well as to show 

that as t — cx), 

P(A;,t) ^ (14) 
in this specific initial state case. 



3 Generating Functions and the Decoherence 
Equation 

3.1 Generating Functions 

The direct calculation involves some formidable, very complicated combi- 
natorics. Therefore, we introduce the idea of generating functions. The 
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generating function of the decoherent quantum random walk is 

oo 

Pm,n{Xi z) — ^ ^ Pm,n{Xi t)z . (l^) 
t=0 

The Fourier transform of the generating function is 

PmAk,z) = J2Pm,n{^,z)e'''\ (16) 



Note that for z in the unit disk {z : \z\ < 1}, since \Pm,nik,t)\ < 1 and 

■it=0 



\Pm,n{^^'t)\ < 1 for every t, the Ylt^oPm,n{k,t)z^'ii and P^,n(a;, 2)'s are ana- 



lytic. Furthermore, 

oo 

$^$^|P™,n(a;,t)e^*^V|2<oo. (17) 



X t=0 

Hence, by Fubini's theorem, we have 



Pm,n{k, Z) = J2 Prn,n{k, t)z', (18) 



t=0 



i.e., Pm,n{k, zYs are analytic and Pm,n{k,tys are the coefficients of in the 
expansions of P^_„(A;, ^)'s. Therefore, instead of finding P^_„(A;, t)'s directly, 
we first find the explicit formulae of Pm,n{k-, ^)'s and then we apply Cauchy's 
Theorem 

for some < r < 1, to get Pm,n{k, t). 

3.2 The Decoherence Equation 

The functions Qni,n{k,z) and Qm.n{x,z) are very important in our proofs. 
We let Wm,n{k,t) = ^^Wjn,n{x,t)e'^''^ be the Fourier transform of the pure 
Hadamard walk. We also let Qm,n{k,z) = ^Yl^i^m,n{k,t){qzY for < 
p < 1 and q — 1 — p. Note that \Wm,n{k,t)\ < 1. Hence, for z G {z : 
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\A < g}) \Qm,n{,k, z)\ < OO. Therefore, Qm,n{k,z) is analytic in {z : \z\ < 
i}. Furthermore, let Qm,n{x,z) = ^'^'^iWm,n{x,t){qzy , and by Fubini's 
theorem again we have 

oo 

Q^,n{k,z) = ^5^ty„,„(A;,t)(gz)* = 5^g™,„(x,z)e^^^ (20) 

^ t=l X 

Using the above notations, we derive the following theorem, whose proof 
could be found in the appendix and in [30] . 

Theorem 3.1 (The Decoherence Equation) 

The functions Pm,n{k,z) 's are analytic in {z : \z\ < 1} and are meromorphic 
in {z : \z\ < Furthermore, if we denote the matrices of {Pm^n{k, z)) and 
{Qm.nik, z)) by P and Q respectively, then 

P=-'il + \l-Q)-\ (21) 
p p 

This equation establishes the relationship between the decoherent quan- 
tum random walk (left-hand side) and the pure quantum random walk (right- 
hand side). By working on the Fourier transform of the pure quantum ran- 
dom walk, we get the formulae of Pm,n{k, z)^s from this equation. 

4 Main Results 

We list our results for the two types of decoherent quantum random walks 
here. The step-by-step mathematical proofs of theorems in this section could 
be found in [30] . 

4.1 Results for the Symmetric Decoherent Quantum 
Random Walk 

The following theorem gives the closed form formula of the Fourier transform 
of the generating function of the walk. This formula synthesizes all the 
information of the walk and is crucial in proving our results. 
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Theorem 4.1 The Fourier transform of the generating function of the sym- 
metric decoherent Hadamard walk, P{k,z), is given by 

P{k,z) 

q{q - cos2 /,•) -2 + p COS /,•,: + (1 - cos h)E (22) 
pq COS kz^ — {pq + p)z'^ + p cos kz + {z^ — 2 cos kz + 1)E 

where 

E = ^{(fz'^ - (1 + cos k)qz + \){q^z^ + (1 + cos }i)qz + 1) (23) 

and the square root in the formula of E is defined through the principal log- 
arithm. 

We first show that the position distribution of the walk is symmetric with 
respect to the origin. 

Theorem 4.2 Let Xf be the symmetric decoherent Hadamard walk on the 
line with < p < 1 and q = 1 — p. Then E{Xt) = 0, Vt. 

We then consider the hmiting distribution. We derive the following the- 
orem for the limiting distribution of the symmetric decoherent Hadamard 
walk. 

Theorem 4.3 Let be the symmetric decoherent Hadamard walk on the 
line with < p < 1 and q — 1 — p. Then the characteristic function P{k, t) 
of Xt satisfies 



P(^,*)=exp(-?i±M±Z^,.)+0(r') 



(24) 



as t ^ oo, i.e., 



^^iV(0/ + ^^^-^) (25) 

in distribution as t oo. 

This theorem states that after a long time, the position distribution in- 
duced by the wave function of the particle is Gaussian. We see from the 
variance of the distribution that it is a mixture of the quantum and classical 
limiting distribution. 

We give a plot of different limiting distributions w.r.t. different p's in 
Figure 1. 



9 




— p=1 

p=0.8 

p=0.1 

p=0.05 

p=0.01 



-20 



-10 



10 



20 



X_t/sqrt(t) 



Figure 1. Limiting position distributions of ^ when p=l, 0.8, 0.1, 0.05 and 0.01. 

We see that when p increases to 1, the distribution is asymptotically stan- 
dard normal, as in a classical walk. On the other hand, when p is decreasing 
to 0, the variance increases to infinity, meaning that ^ does not converge. 

We also find the long-term variance of Xt as follows. 

Theorem 4.4 Let Xt be the symmetric decoherent Hadamard walk on the 
line with < p < 1 and q = 1 — p. Then 



Var{Xt 



p 



2J1 



-t - 



2q^ 



p^/l + q'^ 



P 



- {l + q'-^l + q2) + 0ie 



(26) 



-ct\ 



for some c> 0, as t ^ oo. 



This theorem shows that for fixed p and large t, the magnitude of the 
walk is growing linearly in \/t. 
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Figure 2. Standard deviation of the particle position as a function of p at t=200, 
300, 400 and 500. The grid of p is from to 1 with increment 0.01. 

We plotted the standard deviations obtained from ( l26l) against p's in 
Figure 2. This picture compares very well with Figure 1 in |26j, where the 
values of standard deviations are from numerical simulations. Note that 
when t = 200 and p = 0.01, the standard deviation seems too low. This is 
because when t = 200, the limiting phenomenon has not occurred yet for the 
walk with p = 0.01. Therefore, the formula (1261) is not a good approximation 
and the decoherent walk is still similar as "quantum". We discuss this in 
Section [51 



4.2 Results for the Decoherent Hadamard Walk Start- 
ing at |0, 1). 

Now we consider the decoherent walk starting at |0,1). As before we first 
find its expectation. 

Theorem 4.5 Let Xt be the decoherent Hadamard walk starting at |0, 1) 
with < p < 1 and q = 1 — p. Let fit = E{Xt), then we have fit = 

'^^^'^ ^ -\- 0{e^'^^) for some d > 0, as t ^ oo. 
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This theorem shows that the hmiting expected position of the decoherent 
Hadamard walk is to the right of the origin, if the initial coin state is "right." 
We see that when jo — > 0, //j — ^ cxo. It is consistent with the result in [11] 
that the pure quantum random walk starting with chirality "right" drifts to 
the right. 

For the second moment, we have the same result as for the symmetric 
walk. 

Theorem 4.6 Let Xf be the decoherent Hadamard walk starting at \0,1) with 
< p < 1 and q = 1 ~ p. Then 

Var{Xt) 

P P^^T^ 

2 



-(l + g2-v/rT?) + 0(e-^*) 



P 

for some c> 0, as t oo. 

Now we show that the limiting position distribution of the decoherent 
Hadamard walk starting at |0, 1) is also Gaussian. 

Theorem 4.7 Let Xf be the decoherent Hadamard walk starting at |0, 1) 

with < p < 1 and q = 1 — p. Then the characteristic function P{k, t) of Xt 
satisfies 

>(A.,) = exp(- P+^^-\ ^) + 0(H) (28) 
as t oo, i.e., 



y/t p 
in distribution as t oo. 

Remark 4.1 Note that here the converging speed is 0{t^^) while we have 
0(t^^) for the symmetric walk. This is because when one takes the average 
of the Pm,n{^,t) 's, the error terms in 0(t'2) cancel out one another. This 
result shows that the symmetric walk converges faster. 
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5 Speed of the Walk When p is Small: Pseu- 
doquantum Phenomenon 



In [TT], it is shown that the long-term variance of the Hadamard walk is 
(1 — ■^)^^- We proved the same result via our approach by letting p — >• 0. 

Theorem 5.1 Let Qt denote the Hadamard walk on the line, the long-term 
variance satisfies ^"^2'^*'' ~ as t ^ 00. 

We also investigated the case when p is small. Because current literature 
shows that the distribution of the pure Hadamard walk compares well with 
the uniform distribution over [—-^, ■^], we also compare the decoherent 
Hadamard walk with it. Denote the uniform distribution over [—-^, by 
Ut- We shall compare the variance of Ut and the variance of the symmetric 
walk Xt. Note that by (Ell) and 

Var{Ut) = ^. (30) 
o 

The difference Var{Ut) — Var{Xt) is minimized at 

io=«<V^-^>+3. (31) 

P 

The minimizer p from fl3Tl) when t^ = 200 is about 0.0014. This compares 
well with the numerical result in [26], which demonstrated that when to = 
200, the p that minimizes the difference between Xt and Ut is about 0.0013. 
A plot of to versus p is also given in Figure 3. 
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Figure 3. The pseudoquantum time to as a function of p. The grid of p is from 

0.01 to 1 with increment 0.01. 

The time to is interesting in the sense that the variance of the decoherent 
walk could not be regarded as linear in t before to, i.e, the walk before 
to could not be regarded as "classical." We call the period from to to 
"pseudoquantum" since the walk takes on quantum features. After to, the 
variance of the walk approaches the formula obtained in (l26l) . For example, 
for p = 0.01, to = 247.3. Therefore, when t = 200, the limiting behavior has 
not occurred and the decoherent Hadamard walk has more quantum features 
than classical ones. 

6 Conclusions 

We have investigated the quantum walk with decoherence on both position 
and chirality states. Long-term limits are obtained for both the symmetric 
walk and the walk starting at with chirality "right." We provide analytical 
explanations of the dynamics of the decoherent quantum walk system and we 
see that the system is indeed a mixture of quantum and classical ones. The 
limiting distributions of quantum random walks are shown to be Gaussian if 
decoherence occurs. These results are very important properties of the de- 
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coherent quantum random walks and could be essential for the development 
of quantum algorithms and experiments. 

We also see that when p is small, the system remains non-classical for a 
very long time. If a quantum algorithm can be finished before the classical 
features appear, then we may call it a "pseudoquantum" algorithm. However, 
we do not know how fast the "pseudoquantum" algorithms are as compared 
to the classical ones. Therefore, we suggest future studies on these areas. 
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A Proof of Theorem 3.1 



We start with an observation about the decoherent quantum random walk 
and get a recursive formula. Then we apply that formula to the Pm,n{k, z)^s 
to establish the decoherence equation. 

For any state |0) G H, \(f)) can be written as \(f)) = J2yi{y^^\4>)\yJ)- By 
definition, 

=qP,{\0,m),U*\<P))+pJ2\{yJm'Pmrn),U*\y,l)). (32) 

In particular, for |0) = \x,n), we have 

Pt+i{\0,m),\cj)^,r.)) = Pt{\0,m),U*\x,n)), (33) 
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which in turn gives 

Pt+ii\0,m),m 

=qPt{\0,m),U*m+pJ2\{yJm'Pt+i{\0,m),\yJ)). (34) 

y,i 

This is our recursive formula. Also, for t = 1, we have 

Pi(|O,m),|0)) 

=qmU\0,m)\'+pY,\{y,l\<P)\'P^{\0,m),U*\yJ)y (35) 
y,i 

and 

Pi(|0,m), \x,n)) = \{x,n\U\0,m)\'^. (36) 

Apply the recursive formula (IMIl and (1551) repeatedly, we have the following 
equation 

Pm,n{x,t) = Pt(|0,m), \x,n)) 
= Y.P1^'' E I ^l(^*)'l^' n)\'P,^,{\0, m), \y, (37) 

s=l y,l 

+ g*-i|(x,n|[/*|0,m)|2. 

Note that by the definition of Wm,n{x-,t), we have that 

\{yMUy\x.n)\' 
=\{xMU'\y,l)? (38) 
=Wi^n{x - y, s) 

and that 

\{x,n\U'\0,m)\' = WmA^,t). (39) 
Therefore, (l37j) becomes 

Pm,n {X-i 

= ^pq'~^ Wi,n{x - y, s)Pm,i{y, t - s)+ (40) 

s=l y,l 
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Now, by (HOD, for z e {z : \z\ < 



~^x,n ~^ ^Qm,n{x, z) Qm,n{Xj (41) 

y,i 



where 



c-o.m 



1, a = (3,m = n 
0, otherwise 



Finally, we take the Fourier transform on (HTl) to get 

Pm,n(^k , Z^ 

P I 

where 

1, m = n 



51 



0, otherwise 



The interchanges of summation are justified since the series absolutely con- 
verges. Now, denoting the matrices (Pm,n{k, z)) and {Qm,,n{k, z)) by P and 
Q, we have the following equation: 



which is equivalent to 



P = / + ^Q + Pg, (43) 
V 



P[I-Q) = -^I-Q) + Il (44) 
p p 



We complete the proof by the following lemma. 
Lemma A.l For z & {z : \z\ < 1} , the matrix I — Q is invertible. 
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Proof. For z & {z : \z\ < 1}, if we let Qm,n = Qm,n{k, z), we have 
\Qm,l \ + \Q m,2 1 



<-5^g*(|W^™,i(^,i)l + \W^,2{k,m (45) 



oo 

p 



qp 

5]) implies that ||Q||oo = max^ \Qm,n\ < 1- Therefore, 

oo oo 

II J]Q-''||oo<X1II'5'IIoo<oo, (46) 

3=0 j=0 

i.e., the series Yl^=o converges. This implies that (/ — Q)~^ exists and 

oo 
j=0 

By Lemma lA.ll / — Q is invertible and together with (jH]) we have 

P=-^I+-{I-Q)-\ (47) 
p p 

which is exactly ( l2TI) . 

For z E {z : \z\ < \ det(/ — Q)\ < oo. Hence, det (/ — Q) is analytic. 
Note also that 

A,(fc,.) = -^+ 



Pi,2{k.z) 



p p det {I — Q) 
Ql3 



p p det (/ — Q) 

Therefore, Pm,n{k,zys are meromorphic functions for z E {z : \z\ < ^}. 
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B Proof of Theorem 14. IL 



To get formula (122|) . first we need to know the formulae of Wm,nik,t), i.e., 
we look at the pure quantum walk in the Fourier transform. 

Similar to the setup in [TT], if we let the initial state be |0,m), and 
we let '^rn,n{x,t) = {x , n\U^ \0 , TTi) be the coefficient of the walk at time 
t at coordinate \x,n), then Wm,n{x,t) = |\E'm,„(x, We also introduce 

^m,nik,t) = Ex ^m,n(x,t)e^^" and ^m{k,t) = ( A;, t) , $„,2 (fc, t) )^ in the 
Fourier transform as in [TT]. The evolution operator in Fourier transform 
space, U{k), is defined s.t. 'i/rn{k,t + 1) = U{k)'^rn{k,t). It is obtained in 
[TT] that 

1 / ik ik \ 

m = ^ ( ) . (49) 

Therefore, if we let = i + „ , "^"^ ^ . , and Ck = 777=^^, then for 
t = 2j — 1, we have 

?7*(A;) 

-e-'^'^^ + 2 Ak cos Ukt 2Ck cos uJkt \ (50) 

2Ck cos ujkt e^'^''^ — 2 Ak cos UJkt J' 

Also, for t = 2j, we have 
^*(^) 

g-i'^fct + 2Aa:« sin Wfct 2iCfcsinu;fct \ (51) 

2iCkSmuikt 6*'^*=* — 2^4^^ sin cjfct y 

Now, note that 

^mAKO) = J](x,n|0,m)e*'=^ = (52) 

a: 

and that $m(A;,t) = (f/(A;))*$m(A;, 0), we conclude that $m,n(^,i) = iU\k))n,n 
Hence, for t = 2j — 1, 

^ (k, t)=- e-^^^^ + 2Ak cos Ukt, 
^ i,2{k, t) =2Ck cos Ukt, 

^ _ (53) 

^2,i(fc,i) =2CfcCosa;fct, 

$2,2(fc,t) =e^^'=*-2AfeCosa;fct 
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(54) 



For t^2j, 

$1,1 (A;, t) =6"*'^'=* + 2Aki sin LOkt, 
^lA^jt) =2iCksmukt, 
^2,i{k,t) =2iCk sin Ukt, 
^2,2{k,t) =6^'"'=* - 2AkisincOkt. 

Since Wm,n{x,t) = \^rn,nix,t)\'^, in the Fourier transform, 

^ 1 Z'^'' - 

Wm,n{k, ^) = ^ y *m,n(s, ^)*m,n(A: " S, t)ds. (55) 

We separate the real and imaginary parts of Wm,n{k,t) and get their 
formulae as follows. For t — 2j — 1, 



Re{Wi.i{k,t)) = Re{W2,2ik,t)) (56) 

1 f'^'^ 1 cos LOgt cos LOk^gt /, \, 

= — / cosscos (A; — s)as— 

2% ,/o 2 cos Us cos UJk-s 
r 

sin cOst sin cuk^stds, 



27r 

i?e(l?i,2(A;,t)) = i?e(W2,i(A:,^)) (57) 



27r/ 



cos k cos a;^^ cos cok-gt 



ds, 



27r Jo 2 cos Us cos w^-s 
7m(l?i,i(A;, t)) = -Im{W2,2{k, t)) (58) 



■ coss 



1 r^, 1 

= — / cosa;stsma;jk_si+ 

27ryo 



~i/2cosa;5 

1 cos (k — s) . , , 

H ;= coscjfc-st smtt^stjas, 

V2 cosa;fe_s 

/m(W^i,2(A;, 0) = -ImiW2,iik, t)) (59) 

1 /'^'^ sin k cos uj cos ^'a - s^ 

-as. 



27r Jo 2 cos Us cos c<;jfc_s 
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For t^2j, 

Re{W^,i{k, t)) = Re{W2,2ik, t)) (60) 

1 f'^'' 1 sin Ust sin Uk-st , 

= / cosscos (A; — s)as+ 

27r Jo 2 cos ujs cos oJk-s 



?1-K 



1 r 

H / cosoisi cosa;fe_stas, 

27r Jo 



i?e(W^i,2(A;,t)) = i?e(W^2,i(A;,t)) (61) 
'"^'^ cos sin LOgt sin Wfe-^t 



27r7o 



2 cost's cos a;^-. 



-(is, 



/m(W^i,i(A:, t)) = -7m(W^2,2(A:, t)) (62) 

sin LOgt cos a;fe_st+ 



1 1 coss 



27r Jo cos 



1 cos {k — s) , , 

H ;= sina;fc-st cosa;st)as, 

V2 cosa;fe_s 



/m(W^i,2(A;, t)) = -7m(W^2,i(^, t)) (63) 
"^'^ sin A; sin iOgt sin 



27r7o 



-ds. 



27r Jo 2 cos a;^ cos lo^-. 
Now we are ready to find Pyn^n^i^y^ formulae. We first introduce several 
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short notations. We introduce the T,i^s for z E {z : \z\ < Let 

51 = Re{Qi,,) (64) 

oo 

= ^Y.^Re{W,,,{k,tmqz)\ 

52 = Re{Qi,2) (65) 

00 

53 = /m(Qi,i) (66) 

00 

= ^5][/m(m.i(A:,t))](g.)\ 

54 = /m(gi,2) (67) 



00 



^^J2[Im{W2Ak,tmqzy. 



Since |lVm_„(/i;, t)| < 1, for z G {z : \z\ < i}, the above series all converge. 
Therefore, Sj's are all analytic in {z : \z\ < -}. 
Now det(/ — Q) can be written as 

det (/ - Q) 

= 1 - Ql,l - Q2,2 + Ql,lQ2,2 - Ql,2Q2,l (68) 
= (1-Si)2-S2' + S3'-S4'. 



Note that P{k, z) = ^ J2m n Prn,n{k, z) . By the decoherence equation fl2T]l . 
this function can be written as 

P{k,z) 

p 2p det (/ - Q) (69) 



P p((l-Si)2-E2' + S3'-S4') 

Therefore, once we have the formulae of Sj's, we have the formulae of 
P{k,z). To find Sj's formulae, we first look for the formulae for a real 
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number z E (~qjq)- Then we show that they are the desired formulae for 



all z e {z : \z\ < Let 

oo 

h^Yl - COS [(2j - l)u;k-s]{qz)^'-\ (70) 

oo 

I2=J2 Sin [(2j - l)c^.] sin [{2j - l)uk-s\ {qzf'-\ (71) 

oo 

/3=5^ cos [{2j)us] COS [(2j>fe_3]((Z^)'^ (72) 

oo 

/4=J] Sin [(2j>,] sin [(2j>fe_,](g^)2^ (73) 
i=i 

oo 

/5=^ COS [(2j - l)a;,] sin [(2j - l)a;fe_,](g^)2^-\ (74) 

oo 

76=5^ Sin [(2j - l)ujs] cos [(2i - l)ujk-s]{qzf'-\ (75) 

oo 

/7= J] Sin [(2j>,] cos [(2j>fe_,] (gz)^^ (76) 



oo 



/8=^ COS [(2j>,] Sin [(2j>,_,](gz)2^-. (77) 
Since the Ej's are bounded, we can interchange the integral and the summa- 
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tion to write the Ej's as 
Si 



pi f'^'' .lcosscos(k - s) ^. ^ 



g 27r Jo 2 cos ujg cos Uk 
pi f^"" 1 cos k 



Jo 2 cos a;^ cos ook-s 



g 27r 7o 2 cos Ug cos a;jk_s 
S3 

pi Z"^'" 1 . cos s , ^ COS (A- - .s) , ^ ^ , , , 

q2'K Jo COSCUs COSUk-s 

'■^'^ sin k 



_p 1 r'^ 1 
~q^Jo 2; 



! COS a;^ cos ujks 

Then we have 



-h + /4)(is. 



where 



4 

^-^cosa;^ cos(jA;-sg2;(l - g^2;^), 
-h + h 

: — [— - sin s sin [k — s)qz + g^z^[cos^ s + cos^ (A; — s) — 
3 

— 1] — - sin s sin (A; — s)(f'z^ — q^z"^] , 

h + l7 

- — —z=qz cosa;5[sin {k — s) + 2g2;sin s + q^z'^ sin (A; — s) 
-D \/2 

/6 + /8 

- — —=qz COS ujk^s[sin s + 2qzsin (k — s) + g^2;^sin s], 
D y2 



D 

— cos {k — 2s) {q^z^ — 2 cos kq'^z'^ + qz)+ 
+ q'^z'^ — cos kq^z^ — cos kqz + 1. 
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Therefore, 

Si 

1 f'^'' 1 

—pz — / —[cos {k — 2s)qz{cos k — qz)+ 
2% Jq D 



H — cos k -\ — cos kq^z^ — (fz^]ds, 
2^ -'2. 



pz COS k{l — q z )— J '^^^1 



—pzsm.k— I — [cos {k — 2s) qz + - + -(f z'^]ds^ 
2tt Jq D 2 2 

S4 

1 1 /"^'^ 1 

=-p.sinA:(l-gV)-y^ -rf.. 

By the integral formula 

/dx 2 lb — c ,1 , , 
-; = — , arctan ( \ / tan ( -ax) ) 

6 + ccosax aVb"^ - \b + c ^2 " 

for 6 > c and the fact that 

(l"z^ — cos kcl'z'^ — cos kqz + 1 
>q^z^ - 2 cos kq^z^ + qz 



for 2 e (— i i), we have 



1 1 



1 ^ds 
27r 7o D 



=((1 + qz){l - qz){{q^z^ - {1 + cos k)qz + l)x 
X {q^z^ + (1 + cosA;)g2; + 1))^)"^ 

Letting 

-E = ^ (g2^2 — (1 + cos A;)g2; + l){q'^z^ + (1 + cos + 1), 
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we have 



Si =^-2 ^ -J — [cos k-qz 

q^z^ — 2 cos kqz + 1 



cos k — 2qz + cos kq^z"^ 
2E ^' 

S2=p^cosA;— , (95) 
^ pzsink 1—q'^z'^^ 



q^z"^ — 2 cos kqz + 1 2E 

Now that we have obtained the formulae of Sj's for z G we can 

check easily by taking the principal branch of log, that the formulae are 
analytic in {z : \z\ < ^}. Hence, by the Analytic Continuation Theorem, 
they are the desired formulae for2;G{2;:|2;|<^}. 

Finally, the theorem is obtained by applying fl95|) to fl69l) . 



C Other Proofs in Subsection 14.1 



Proof of Theorem \4.2\ Note that from the formula, for some r < 1, we have 

EiXt) 

1 / ^^^(0'^)^^ ^^^^ 



2m J\z\=r iz^~^^ 
=0. 

The change of the order of integration and differentiation is justified since 
dkP{k, z) is continuous on the contour. 

Proof of Theorem \4.3\ The denominator of P{k, z) has less than eight isolated 
roots. We shall now look for the root with the smallest absolute value. This 
root has no closed form. However, since we concentrate on the asymptotic 
behaviors, we need only to know its behaviors around A; = 0. The properties 
of this root are summarized in the following lemma. 
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Lemma C.l Let D{k,z) denote the denominator of P{k,z). Then the root 
of D{k, z) = in z, with z = 1 when k = is of the smallest absolute value 
in a neighborhood of k = 0. If we denote it by z{k), then z{k) has multiplicity 
one and can be written as follows. 

z{k) = l + dkz{0)k + o{k). (97) 

Proof For k = 0, D{0,z) = (1 - z){l - qz){pz + (1 - z)^/lT^). By 
solving this equation we can see that z = 1 has the smallest absolute value. 
The root of the second smallest absolute value has a closed form expression, 
which can be found in the appendix of [3D]. We denote this root by z{p). An 
expansion of the root around p = is 

~z{p) = 1 + + 1(1 + -Ly + o{p^). (98) 

Now, by continuity of k, z{k) has the smallest absolute value in a neigh- 
borhood of = 0. 

Since d^D^k, 2)1^=0,2=1 7^ 0, z{k) has multiplicity one. We then apply the 
Implicit Function Theorem to find its derivatives. 

Remark C.l For p 1, D{0,z) ^ 1 — z, for all z E {z : \z\ < which 
implies that other roots go to infinity. In the limit, when p = 1, there is only 
a single root. 

Now we utilize the Implicit Function Theorem to find dkz(0) and dk^z(0) 
where z{k) is defined implicitly by D{k,z) = 0. 
Taking the first derivative we get 

0=dkDik,z) 
= — pq sin kz"^ + pq cos k3z'^dkZ — (pq + p)2zdkZ— 

— p sin kz + p cos kdkZ+ (99) 
+ E{2zdkZ + 2 sin A;^ — 2 cos kdkz) + 
+ dkE{z^ -2coskz + l). 

For A; = and z = 1, the equation becomes 

{pq - p)dkZ = 0, (100) 
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which imphes that 

dkz{0) - 0. 

Now, for dk^z{0), we can take the second derivative on D{k, z) to 

— — pq cos kz"^ — 2pq sin k3z'^dkZ+ 
+ pq cos k{3z'^dk'^z + 6z{dkzy) — 

— {pq + p){2zdk' z + 2{dkzY) — pcoskz— 

— 2p sin kdkZ + p cos kdk^z 

+ E{2zdk^z + 2 cos kz -2 cos kOk"^) 
+ 2dkE{2zdkZ + 2 sin — 2 cos /cSfez) 
afe2£;(^2 -2cosA;z + 1), 

which in turn gives 

Similarly, taking the third derivative of D{k, z) — Q gives 

dk^z{{)) = 0. 

Also, by taking the fourth derivative we get 

9,S(0) ^^_^(76g^ - 83g^(l + 

+ 16g-'^ + 68g2 - ^^(1 + q^f^ - 37g(l + + 
+ 16g-23(l + g2)5 + 28). 
Hence we have the expansion of z{k) at /c = 



,(,) = i + £±M±Z^,^ + o(*^). 

2p 

The residue of ^^^^ is 

We now prove another lemma. 
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Lemma C.2 



lim {z{k) - z)P{k, z) = l + 0{k^) 



as k 0. 

Proof. Note that 7^ 1, 



i.e., Ve > 0, 36, s.t., 



hm{z{k) -z)P{k,z) = 1, 



\{z{k) - z)P{k,z) -1\ < e 
for I A; I < S. flllOp implies that 



for I A; I < 6. Hence, 



for |A;| < 6, i.e.. 



hm \{z{k) - z)P{k,z) - 1| < e 

2— >2(fc) 



hm (^(A;) - z)P{k,z) - 1| < e 

Z— >2(fc) 



hm hm (z(k) — z)P(k, z) = 1. 

fc^0 2^2(fc) 



Now, for a smah ri > s.t. z{k) is the only pole inside the circle 
we have 

hm ]im Uiz{k)-z)Pik,z)-l) 

k-*Oz^z{k) k 



\mv\{— i P(k,z)dz-l) 



2ni 
-0. 



dkPiO,z)dz 

\z-l\=ri 
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Similarly, we also have 



lim lim ^((z(k) - z)P(k, z) - 1) 
f dk'^P{0, z)dz 

J\z-l\=ri 



~^^^^{l-zY {l-zY{pz+{l-z).JliqH^y 

Therefore, \Ym^^^^k){z{k) - z)P{k, z) ^ 1 + 0{k'^). 

Now for any fixed k e [0, 27r], the characteristic function of ^ is P(;^, i). 



< 



Since the roots of D(k, z) are isolated, we can set r(p) — 1 + s.t. \z(^ 
r{p) and other roots are outside the circle {\z\ = r{p)}. Furthermore, when t 
is large, is small, hence the lemmas are applicable. We define the contour 
C as C = {z : \z\ = Tq} [j{z : \z\ — r(p)}, where Tq < 1. 



By definition. 



P(^,t) = _L/ ^^dz. (116) 



\z\=ro 



Z 



Since is the only pole inside the contour, we have 

-«e.(%l,4), 



(117) 



'|z|=ro ^ ''|z|=y(p) 

For fixed < p < 1, supfe^|^|=^(p) |P(^, z)| < oo. Hence, 



|z|=r(p) Z 



dz^O{r{py'). (118) 



We have 



P( A, i) ^ _i?e5(^^, ^( A)) + 0(r(p)-*). (119) 
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Note that by ffTUg]) . we have 



hm hm (2(-^)-z)P(-^,2) = l + 0(ri). (120) 



Note also that by ffT06|l . 



.(i.)^l, ^^^^-^ g,0(n. (121) 



which imphes that 

, 1 



\t+i 



p+2v/rT?-2g 

2p 



(122) 



Therefore, by ffTOTD . V/c G [0,2 



TT 



P(.( A), = exp{-^A^4^^P} + O(r^), (123) 

as t oo. Hence, the hmiting distribution of the symmetric decoherent 

Hadamard walk is Gaussian with variance v = P+^V^+g — 

p 



Proof of Theorem \4.4\ For Xt, we can also find its long-term variance. Let 
C be the same contour as before, when t is large, 2; = 1 is the closest root to 
among all that of the denominator of P{k, z). 
Note that 

- dlP{Q, z) 



^ 2z^{-l + v/TTgV) (124) 
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and that 



2(v/rT^-l), , 2g2 2 



(125) 



Hence, 



-1 - ^^^^^^^ + — =^ + ^(1 + - v^TT? 



= -afc2p(0,t) 

1 / dk'P{0,z) 



2Tci Tc z^+^ 



dz 



{I - zfipz + {I ~ z)^l + qH^y 
p + 2y/l + q^ -2^ 2g2 



(126) 



P PVl + q' 

- 4(1 + - + g2) + 0(r(p)" 



The change of the order of integration and differentiation is justified since 
dk^P^k, z) is continuous on the contour. Hence, the long-term variance of 
the walk is WpZl^^ _ _2gi^ _ + ^2 _ yiT^) + 0(r(p)-*). 



D Proofs in Subsection 14.2 

From the decoherence equation we have 
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Let Xt be the walk starting at position with coin state 1 at time t and 
P(A;, z) = Pi iik, z) + Pi, 2(^5 z) be its generating function. Then 

P{k,z) 

_ q 1 - Si + iSa + S2 + iT^i 

~~p^ pdet (I-Q) (129) 

='^^^'^^ + Vdet(/-g) - 

Note that Sc, and S4 both have a factor of sin k, we denote and 

_ _ ^ ' sm k sin /c 

by S3 and S4 respectively. 
Proof of Theorem \4.5[ Note that 

., S3 + S4 „ (130) 

=H )lfc=o 

P 



(1 - z){pz + (1 - 2) v^TTg^^) 

Let C be the same contour as before. When t is large, z = 1 is closest to 
among the roots of the above denominator. Hence, for fixed p, 

E{Xt) 

^-9fcP(0,t) 



1 r dkPio,z] 



^ dz 
27ri J(j iz^~^^ 

± / iVlTW'-l) (131) 



27ri Jc - z){pz + (1 - z) 

2*(1 - ^)(pz + (1 - 2) yiTg^^) 
=^^^^^ + 0(r(p)-). 
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Proof of Theorem \4.(\ Note that — (9^P(0,2;) must be real. Hence, 

-dfp{Q,z) = -dlP{^,z). (132) 
Therefore, the formula is the same as before. 

Proof of Theorem |^. 7[ Let Ht = E{Xt). We want to show that 
N{0,v). The long-term characteristic function is 

^, k . _iiHh. 

, , , , V . (133) 

-iiiii . fc 1 / 1 , S3 + E4 ., , , ^ ' 

Hence, the limiting distribution of the decoherent Hadamard walk start- 
ing from with coin state 1 is Gaussian with variance v = 

o p 

E An Alternative Proof of Theorem 15.11 : 



By Theorem 14. ![ 



-dlP{0,z) 

z 2z2(-i + ^1 + g2^2) (134) 



(1 - (1 - z)\vz + (1 - 2) v/TT^) ■ 

Note that the variance of the walk at time t is given by 

Var{Xt) = -dlP{{),t) (135) 

and 

P(0,z) = J]F(0,t)^*. (136) 
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Thus, Var{Xt) is the t-th coefficient of the Taylor expansion of — 9^P(0, z). 
Asp ^ 0, Xf becomes Qt- In particular, the function Vo{z) = J2t=i Var{Qt)z^ 
can be obtained from I134[ 

{i-zy [i-zy vi + z^ 

Comparing the coefficients of the Taylor expansion of Vq{z) gives 

VariQ,) = ^ - E - 2j)(t - 2j - (138) 



As t — *• oo, 



VariQt 



•2^ n(j!) 
1 
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